February 7, 2008 14:52 WSPC/INSTRUCTION FILE OptionPriccCou- 
pledPriceTime 



repetowp&tcd.ie 
mcubed@unr. edu 
richmond(§tcd. ie 



February 7, 2008 14:52 WSPC/INSTRUCTION FILE OptionPriccCou- 
pledPriceTime 



International Journal of Theoretical and Applied Finance 
© World Scientific Publishing Company 



Pricing of options on stocks that are driven by multi- dimensional 
coupled price-temporal infinitely divisible fluctuations. 

Przemyslaw Repetowicz 
Department of Physics, Trinity College Dublin 2, Ireland 

Mark M. Meerschaert 

Department of Physics 211 Leifson Physics Building University of Nevada Reno NV 8955 

Peter Richmond 
Department of Physics, Trinity College Dublin 2, Ireland 

Received (29/09/2004) 
Revised (Day Month Year) 

We model the price of a stock via a Langevin equation with multi- dimensional fluctua- 
tions coupled both in the price in time. We generalise previous models in that we assume 
that the fluctuations conditioned on the time step are assumed to be compound Poisson 
processes with operator stable jump intensities. We derive exact relations for Fourier 
transforms of the jump intensity in case of different scaling indices E_ of the process. 
We express the Fourier transform of the joint probability density of the process to attain 
given values at several different times and to attain a given maximal value in a given 
time period through Fourier transforms of the jump intensity. Then we consider a port- 
folio composed of stocks and of options on stocks and we derive the Fourier transform 
of a random variable S)t (deviation of the portfolio) that is defined as a small temporal 
change of the portfolio diminished by the the compound interest earned. We derive a 
functional equation for the price of the option on the stock as a function of the stock 
subject to the condition that the deviation of the portfolio has a zero mean E [Si] = 
Therefore the option pricing problem may have a solution. 

1. Introduction 

Early statistical models of financial markets assumed a Gaussian distribution J^. 
However, there is evidence 0? that stock price returns are not Gauss distributed and 
that the distributions diminish slowly in the high end (fat tails). A better statistical 
description is provided by a model where the logarithm of the stock price is a one 
dimensional compound Poisson process with jumps being asymptotically power-law 
distributed (stable Levy process ^SJ). That model, however, does not eliminate dis- 
crepancies between experimental facts and the theory and drawbacks of the theory. 

The exponents of the power laws ( tail indices ) in distributions appear to be out- 
side of the range of validity predicted by the theory '1 1. Furthermore models in the 
literature assume that jumps in the price are independent of waiting times 
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between jumps (no price-temporal coupling). 

Fat tails can be explained by models basing on operator stable Levy processes 
J^. Here the price of the stock is an exponential from a sum of a large number 
of projections of many-dimensional jumps whose distributions are invariant under 
auto- convolution. This implies that the distributions are, in the high end, mixtures 
of power laws with different exponents, each one describing a walk in one of the 
dimensions. 

We briefly describe development of models that account for a price-temporal cou- 
pling ( coupled random walks ) below. 16^7^ formulate an analytical model that yields 
probability distributions of price jumps under an "ad hoc " ansatz about the form of 
price-temporal coupling. '8! derives a long time scaling limit for the jump probabil- 
ity distribution in coupled random walks for an anomalous fractional diffusion. ^10! 
obtains results for master equations of coupled random walks in the long time scal- 
ing limit and show their connection to fractional derivative calculus. 11 If considers 
one- dimensional continuous-time coupled random walks (Weierstrass flights) and 
derives the long time scaling limit of the mean squared displacement of the walker 
as a function of two parameters a and (3 which are fractional dimensions that repre- 
sent the flight and the waiting. An extension of Weierstrass flights to a model where 
the walker, after random waiting times, moves with arbitrary finite velocities to new 
positions is formulated in and formulas are derived for the sojourn probability 
density and for likelihoods of extreme events. Finally f ^516^ construct contingency 
tables between logarithmic price returns and waiting times for stocks on the New 
York Stock Exchange and for United States Dollar - Japanese Yen exchange rates 
and reject the hypothesis of the independence of returns on waiting times at a high 
confidence level. 

We formulate a multi- dimensional random walk model that recognises that the 
jumps are not independent of the waiting times. The model has a two-level struc- 
ture with the higher level being a walk with random fluctuations preceded by random 
waiting times and coupled to those and the lower level describing the fluctuations 
as sums of random numbers of random jumps with a jump intensity (Levy density) 
vj{x\t). The levels of the model have a subordinated time horizon that will be ex- 
plained in the next section. Our model works with two arbitrary functions, viz the 
Levy density w{x\t) with a scaling index E_ and a waiting time probability density 
function (pdf) pAxit) with a scaling index < /3 < 1 and generalises previous 
models proposed in the literature in this subject matter. 

2. The model 

In this section we formulate an assumption about the price evolution of the stock 
and we recollect some definitions from the theory of stochastic processes. In section 
\2.1\ we recapitulate some known facts concerned with limit distributions of sums of 
independent random vectors In sections \2. ^ and \2. !A we obtain exact relations 
for Fourier transforms of jump intensities and for Fourier transforms of distribu- 
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Hons of integer powers of one dimensional projections of jumps that are sums of 
large numbers of independent random vectors. In section \2.4\ we derive the Fourier 
transform of a joint probability density of an event that the stochastic process attains 
given values at given times and, in addition to that, that it attains a given maximal 
value in some predefined time interval. In section \2. ,5)) we define an option on the 
stock, define a portfolio that is composed of stocks and of options on stocks, define a 
random variable ( deviation of portfolio ) that describes the temporal change of that 
portfolio in a small time span diminished by a compound rate of interest earned on 
the portfolio in that time span, and we derive conditions on the deviation of the 
portfolio to have a zero mean. 



2.1. The stock price 

Let log(S't) be the logarithm of the price of the stock (log-price) at time t. We assume 
that the log-price has a drift a. In a small time interval St the log-price may or may 
not change with respect to this drift. The probability of that change depends on a 
presumed waiting time pdf pAxit). This means that: 

I with probability 1 — Jq PAT(^)d^ 

where P is a real number that will be specified later, the parameters a and a :~ 
(cTi, . . . , (Td) (volatility) are assumed to be constant as a function of time t and of 
the stock price St. We assume that the random vector 6Lt (fluctuation) is infinitely 
divisible (see e.g. \3l for the precise definition). This implies that the logarithm of 
the Fourier transform (log- characteristic function) of it has a unique representation 
1^ in terms of the Levy-K hint chin formula. This means that: 



X,iS^) := E ^exp(zfe • bU)\ = ia-k--Q{k,t)+ (^e^^-^ - 1 - Yljq^ j 0(f |t)df (2.2) 

where a is a constant, Q{k,t) is a quadratic form that describes Gaussian fluctua- 
tions and the integral represents a compound Poisson process with some rate c > 
and with a jump pdf (j)(x\t). Here we do not use the measure-theory notation but 
instead we replace the measure by Riemann integrals with some probability densi- 
ties. We denote SLt = IT> [a, Q, (j)]. Since in our model \2.1]) the drift of the stock 
price is already accounted for by the parameter a and we neglect the Gaussian part 
Q{k,t) in the fluctuations. We also assume that (/)(x|t) is symmetric in x. Therefore 
we have: 

SLt = IV [0, 0, (f)] (2.3) 

The probability density function of fluctuations conditioned on the random waiting 
time ST reads: 

P [6Lt = Sh\5T = Stj = ojgi^^^^ (^5lt\St'^ = exp(-c) ^ ^(l)''®{6l't\St) (2.4) 
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[^J The joint pdt ot a fluctuation 
SLt — Sit preceded by a waiting 
time 6T — St (the joint fluctua- 
tion pdf) 




[B] The pdf of a fluctuation SLt — 
Sit conditioned on a random wait- 
ing time ST = St (the conditional 
fluctuation pdf) 


^SLtlsri^^tlSt) 


[C] The joint pdf of a jump S£,t — 
Sit preceded by a random waiting 

Llllie U-i Uty 1 LllC JOlllL J uiiip ptii ) 


HSk,St) 


[D] The pdf of a jump S£,t condi- 
tioned on a random waiting time 
ST = St (the conditional jump 
pdf) 


HSklSt) 


[£] The pdf of the n*^^ power 
S(") (^a-SLtY - ^ of the 
fluctuation term in H2.1|l condi- 
tioned on the waiting time ST — 
6t (the conditional n-marginal 
pdf) 


i/(")(z|<5i) 


[J-'] The Levy measure of jumps 
(the jump intensity) 


w{Slt) 


[Q\ The pdf of a waiting time 
ST = St (the waiting time pdf) 


PsriSt) 



Tabic 1. Summary of definitions 



where rt® means an n-times auto-convolution. The last equality in represents 
a probability density function of a compound Poisson process with rate c and jump 
intensity (j){Slt\St). 

The joint pdf ujgi^ grp{Slt,St) of a fluctuation SLt = Sit preceded by a waiting 
ST = St reads: 

^SLuSri^^t, St) = ujgi^^gr^ {sit\St^ PsriSt) (2.5) 

Before proceeding further we summarise the definitions in Table 
Note that the quantities [A], [B] differ from their counterparts [C], [D] in that, 
since the fluctuations are compound Poisson random variables, a fluctuation SLt 

consists, in general, of several jumps 5Lt ~ Ss!f'^ where the number of jumps 

Nt is a Poisson random number Nt = Poisson(c) with rate c. Therefore the quan- 
tities [A], [B] are expressed via sums of auto-convolutions of [C], [D] according to 
formula \2.4\j - We also stress that all statements in this paper are concerned with 
large x scaling limits (the high-end of the distribution) of distributions in question. 
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We assume that the time horizon ST and the price fluctuations SLt are small 
at the higher level of the model and large at the lower level. Note that such an as- 
sumption, accompanied by additional ad hoc hypotheses, was tacit in previous works 
(equation (16) in ^ cond-mat/0308017 preprint). This assumption will enable us 
to use the scaling limit of the jump probability density function in a coupled random 
walk: This is the Theorem 2.2 on page 733 in '81 that reads: 

Theorem 2.1. Assume that the waiting time is (3-stable with index < /? < 1, 
meaning that the Laplace transform of the waiting time pdf has a following form: 
Ct [p] (s) = exp (-/Cr(l - (3)sP) where K. = IC{a) is a normalisation constant of the 
waiting time pdf (JC{a) — p{t)dt) for some a >> 0. 

The large x and large t scaling limit of the of the joint jump pdf in a coupled 
random walk reads 

cj>{x,t) = det{t-''&)w{t-^Ms)'^^ (2.6) 

where vu : R'* B x ^ 'cu{x) G K+ is the probability density in x (termed as the jump 
intensity), and the matrix E_ is an exponent of the operator stable law (stable index) 
corresponding to the jump x G R'^ random variable and det(t^^=) is a normalisation 
constant of the jump pdf. 

The eigenvalues X of E_ satisfy Re[X\ > 1/2 and the eigenvalues Re[X\ = 1/2 are 
simple (not degenerate). 

Corollary 2.1. The scaling limit of the conditional jump pdf is: 

^("^1^) = f ^If^'L- = det{t-P^u.{rmx) (2.7) 
J (p[x,t)dx 

This follows in a straightforward manner by integrating over x in 1^2. b\) . 

Corollary 2.2. The jump intensity zu satisfies a functional equation: 

^E(^Um^^ ^ tn*^(x) = w{t-M=x)det{t-=) (2.8) 
= MO 

where ■uo*'® is the tth auto- convolution defined as a limit of t„-times auto- 
convolutions for some rational tn such that limn^ootn — t. Since the random process 
is infinitely divisible the definition is correct. 

The corollary follows from the self- similarity of the random walk. This means 
that there exists a sequence of linear operators Bin) G L(R'^) and a sequence of real 
numbers 6„ and a sequence of integer numbers kn such that for some independent, 
identically distributed random variables Yi G R'^ we have: 



Yi + ...Yk^] - bk„ Y with a pdf vo (2.9) 
On the other hand if we take the sequence of rational numbers tn we have 

B{n) f fi + . . . ft„fe„) - 6t„fc„ ^ fi with a pdf ro*® (2.10) 
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Now if we take n large enough we can write fc„i„ = fcm(n) for a certain subsequence 
m{n). Therefore the random variable Yi has the same pdfzu except that the argument 
of the pdf is linearly transformed x — > t~=x. 

Corollary 2.3. The Fourier transform ^^{k) := !Fg[w\{k) of the jump intensity 
w satisfies a functional equation: 

w\k) = witM^k) (2.11) 

Where is the transpose of E. The right hand side of the equality 1^2.11]) follows 
in a straightforward manner by taking a Fourier transform of the right hand side of 
\2.f^) . The left hand side is also straightforward. 

Ts [w^®] (fc) = limn^ooJ'x [cc7*"®] (fc) = Zim„_,ooW*"(fc) = ^\k) (2.12) 

see for details. 

We illustrate Theorem 2. 1 in following examples: 
Example 1. Let d — 1 and Yt = Yt — Normal{0,t) is normal with mean zero and 
variance t. Assume that P =. ThenE_— E~ 1/2 andm{x) = (27r)^"'^/^ exp(— a;^/2). 
Example 2. Let d = \ and let the log- characteristic function of Yt be 
log(-E[exp(«fcy)] = — t|fc|^ where < /i < 2. Then E = l/ii and the function w 
can has an asymptotic expansion: 

, N 1 . ,7r ,r(/in + l) 1 ,„ , „x 

n=l 

Ln the next section we are going to analyse the Langevin equation 1^2.1]) for the 
log-price of the stock driven by operator stable fluctuations. For this purpose we 
need the time dependence of the conditional jump pdf on the waiting time for small 
waiting times and we need the probability distribution of the n^^ power of a scalar 
product S'"' :— (a ■ 6Lt ) for n £ N (the conditional n-marginal pdf). We have: 

Proposition 2.1. The conditional jump pdf (f>(x\t) depends on time t as follows: 

(jj(x\t) = t'^^''=zu{t-'^^x) (2.14) 

This is a straightforward conclusion from \2.1\j and from the fact that 
det 



Tr\A\ 

- for every linear operator A G L(IR ) 



Proposition 2.2. The Fourier transform of the conditional n-marginal probability 
density function v'^'^\k\5t) J-z [j^^"^] (fc) as a function of the Fourier transform 
of the jump intensity ■uj{\) := T3 \^\ (A) takes following form: 



oc 

J w\ 



(2.15) 
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where 



/C'"'(t,[) := 



2Re 



2e*^ (j^e~^^" cos |^cos(7r ^"'2,^ e?;en 
and 6 := fccr". For small values of k the conditional n-marginal pdf reads: 



if n is odd 



where 



C5 := 



2cos(.^) 



(2.16) 



(2.17) 



(2.18) 



and £ is defined in Appendix A The proof of the proposition is in Appendix A. 

Corollary 2.4- The mean value of the n^^ power of the fluctuation term in \2.1\) 
conditioned on the waiting time ST = 6t is infinite if n is even and n > 2 and is 
zero if n is odd. 



E 



SL, 



ST = St 



00 for n G 2N and n > 2 
for ne N\2N 



(2.19) 



The first statement follows from diy^'^\k\St) / dk\ = 00 (see \2.1'J\j ) and the sec- 
ond statement follows from the assumption that the conditional jump pdf (l){x\St) 
and hence the jump intensity 'uj{x) are symmetric in x (see beginning of section 
flTTl) . 



Corollary 2.5. The n^^ power of the fluctuation term in 1^2.1]) has a fo 



[a-SLtY ^y\{Stf + 0{{Stff^) 



(2.20) 



where d\ is a random number that depends neither on time t nor on St. 
We have: 
1 

dz 



-P [z < Z < z + dz \ ST = St) = v'^''\z\St) = 



(^) 



(2.21) 
(2.22) 
(2.23) 
(2.24) 



[1 - (Stfi^ - 1) + O {{St)^")] (z) = 

S{z - zo) ~ {StfT^^ \{m - 1) (©fci/"^)] (z) + O {{Stf^) 

for some zq G M. \2.24^ is equivalent to \2.2(]\) for z > zq. The random variable £H 
conforms to a distribution whose characteristic function is (to — 1) (©fc"'^/"(7) . 

The concrete dependence of v^'^'^ {k\St) on time St depends on the form of the 
matrix E_ (stable index). We investigate that dependence in next sections. 

Now we are going to investigate the functional form of the jump intensity 'cu{x) 
in the limit of large jumps x. It is more convenient to analyse the Fourier transform 
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uj{k) for small values of k. This function is a solution of equation \2.11]) and is 
determined in a unique way via the stable index E_ G L(M'*). Since, however, on 
the grounds of the Jordan decomposition theorem, every matrix E_ has a unique 
representation as a block diagonal matrix where every block is of the form: 



(2.25) 











... 


\ 












... 


\ 


1 


a 





... 






/ 







... 







1 


a 






or 













Vo 






1 a 


) 




I 






/ » 


/ 



where a is a real eigenvalue of E_ in the first case, and in the second case 



53 



a -b 
b a 



and I ■■ 



1 
1 



(2.26) 



where a ± ib is a complex conjugate pair of eigenvalues of E_ the set of all possible 
jump intensities w is narrowed down to few classes of solutions only, each one 
corresponding to a particular Jordan decomposition of the matrix E_. In the following 
we investigate these classes of solutions as a function of E_. We proceed as follows: 

(1) Take a Jordan decomposition of E_ G L(M'*), 

(2) Find the Fourier transform w(k) of the jump intensity by solving the functional 
equation Y2.11]) 

(3) Find, from 1^2. 15\) the Fourier transform (A:|(5t) of the conditional n-marginal 
pdf (see Table^^ as a function St and of k for small St and small k. 

In the following sections we seek for the solutions vj of the equation \2.f^ for dif- 
ferent types of blocks in the Jordan representation. The names of the sections refer 
to the type of Jordan decompositions of E. 



2.2. The pure scaling 

LetE= [dfiy^I wh ere I is ad dimensional identity matrix and fi > is a constant. 
Here TrE_ = . The mapping: 



t= -.W^ 3 k^t^'^'^'^ ^kcW^ 



(2.27) 



changes the length of k only. We therefore assume the jump intensity w(k) = Tu{\k\) 
to depend on the length of the vector only. Substituting that assumption into Y2.11]) 
and taking t^'^'^) — \k\~^ we get: 

w{\k\f\''"''' 



and 



w(k) — 



^w{\k\-^\k\) 
exp (c\k\'^^'^ 



(2.28) 
(2.29) 
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where C :— log(ti7(l)). 

The conditional n-marginal pdf (see Table^ from ^2.1 5}} reads: 

POO 

v^''\k\5t)^ d[exp{C(Jt)^|I|''^}/C(")(fc(T",l) (2.30) 



where the kernel JC is defined in \4-ll\l and has a series expansion in I given in 
{4-l^ - From {2.17)) we get the n-marginal pdf for n even and for small values of k. 
We have: 

^(")(fc|5i) = expjc ( pLyiSifk"^] 



V 2an / 

where a„ := cos(7r-^^^^) and the constant £ is defined in Appendix A. 



2.3. The scaling & rotation 

Take d — 2 and let: 

/(2Ai)-i -b 



(2.31) 



I-(^7 .o::-i) (2.32) 



The trace TAe] = n ^. We denote by O ^ :— ( '^^[^2 ^' [^2 I a two dimensional 

y=l3y sin(/?) cos(/3) J 

rotation by an angle (3. The mapping: 

tl"^ : R2 9 fc ^ , , ,,fceR2 (2.33) 

=-blog(t) ^ ' 

changes the length of k by a factor t^'^'^^ rotates by an angle —h\og{t). 
We seek solutions of \2.11]) in the form: 

6j{k) =eyi^[y{\k\,6)^ (2.34) 

where the function v{k) = v{\k\,9) is assumed to depend on the length \k\ and on 
the angle 9 := <(fc, (1,0)). Inserting {2.34^ into {2.11}) we obtain: 

ti{\k\,e) = v{t^^f^''\k\,e - b\og{t)) (2.35) 

Now we assume that the function v factorizes: 

v{k) ^ di{\k\)Q{e) (2.36) 

Inserting \2.ilb]) into \2.il5\) and taking t'-^'^^ ^ = |fc|™ with m eN we get: 

l^p^"ui(|fc|)e(6i) = vi{\kr\k\)eie-b\og{t)) = vi{\k\"'+'^)e{e^2fibm\og{k))i2.37) 

After separating variables and assuming that the jump intensity 'cu(x) is even in x 
we obtain v in the following form: 



}(fc) = ui(l)fc2M Re 



exp 



/i5 niq log(fc) 



(2.38) 
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where {nq,mq) G N for g G Z are some numbers that depend on t. We find only 

solutions for t — |fcp'^"' for to G N and leave the generic case for future work. 

We check that is indeed a solution of 1^2. 'J5\) for t = jfep^™ under following 
assumptions: 



V 



Ipi 



= 2N 



q>l p'GN TOg(TO+l) mpi TOgTO+1 

This leads to a following solution for the jump intensity in the Fourier domain 

/27r (m + e 



(2.39) 



w(fc) = exp 



Vk^^' cos 



\ /i6 m log(fc) 



g— — oo 

where T) := ui(l) For small k the conditional n-marginal pdf reads: 



i>(")(fc|<5t) =exp 
where 



2pi 



P(©crfcl/"j {5tf COS 



q— — oo 



27r (to + 



n{5t,a, k) 



e~bp\og{st) 



n{6t,a, k) 



log(G5CTfci/") + |-log(,5t)^ 



2.40) 



(2.41) 



(2.42) 



where the number © is defined in \2.1t^) . We summarise the results for the Fourier 
transforms of the jump intensity and of the conditional n-marginal pdf in both cases 
of pure scaling and of scaling & rotation: 



vj{k) ~ 



i>'^''\k\St) = 



expfcifcl''^ 

Pfc2, g eos^^il2±ii^ 



cxp 

exp 
cxp 



log(fe) 



pure scaling 
scaling & rotation 



{C(©afci/")''^ (Stf} 



(2.43) 



pure scaling 

(2.44) 

scaling & rotation 



where (5 is defined in ^2.11^) and the results for the conditional n-marginals hold for 
small k only. 



2.4. The joint transition probability 

The purpose of this section is to derive probability distributions of financial instru- 
ments like barrier European style options with multiple exercise times and American 
style options (exotic options) '15] as a function of of the conditional I -marginal pdf. 
We derive a joint probability distribution of log-prices Xt — log(S't) of the stock at 
certain times and of the maximal log-price in some time interval. 
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We fix I + 1 values of time to < ti < t2 < ■ ■ ■ < ti and we analyse the joint 

cumulative distribution function Fx^^.Xt^ Xt, ;A/t, (a;i, a;2, . . . , a;;; y) related to the 

log-prices Xt- for i = 1^ ... ,1 of the .stock and to the maximal value of the log-price 
Mt := SMPo<s<(Xs. We set Xtg = 0, we define St = {ti — to)/N and we define 
integers 1 < ii < . . . < ii < N such that tj — <o = {5t)ij . From the definition of the 
joint cumulative distribution function and from the the SDE \2.1\) we have: 

Fxt^.Xt^,...^Xt,-Mt, (xi,X2,...,xi;y) := P {Xt, <xi,...,Xt, < xf, Mt^ < y) (2.45) 

= p I y Xt,+,^^st) < a^,; U ^to+jm < y] (2.46) 
\j=l J=l / 

= U X! ■ ^^^to+pst < (xj - tjn); IJ X! ^ ' ^^^to+pst j < (y - (io + jSt)n) 
\p=i / j=i \p=i / 

(2.47) 

N 

= J lA(.i,...,.,;,)nKejl'5^)^^J- (2-48) 

I 

where [J denotes a union of I conditions, the integration region 

3 = 1 



A(xi,...,a:,;2/) := <( V (2.49) 

j=l 

is an infinite region in bounded by I -\- N hyperplanes of dimension N — 1 and 
^(^10 ■= '^-•^ll,^^!^) conditional probability density of the random variable 

a-SL. The Fourier transform of that probability density is given in \2.44\l - Now the 
procedure consists of following three steps. 

(a) Derive the joint probability density function fxt,,Xt^,...,Xt^:Mt^{xiTX2T ... ,xi;y) 
by differentiating the cumulative density function, 

(b) Derive the joint characteristic function by taking Fourier transforms of all 
I + 1 random variables and expressing the result through Fourier transforms 
v^^\k\5t) := Tz \v'^^^{z\8t^ (k) of the pdf of the random variable a ■ SL 

(c) Derive the joint pdf in (a) by inverting expressions containing the Fourier trans- 
forms iy^^^klSt) in obtained (b). 

We describe the steps in the following and give examples. 

(a) 

fxt,,Xt^,...,Xt^;Mt^{xi,X2, . . . ,xi;y) := (2.50) 
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Ql+1 r 1 

dx^...dx,dy r^*..^*..-A,;M,(a^i,^2,...,x,;y)J (2.51) 

IN./ \ / i \ N 

j=i \ p=i y V p=i / 

The sum on the right-hand side in \2.5}!^) corresponds to intersections of all 
possible pairs of hyperplanes that bound the region \2.4y^ - 
(b ) The joint characteristic function of the I + 1 random variables in question reads: 

XXt^,Xt^,...^Xt,-.Mt^{ki,k2, . . . ,ki;w) (2.53) 
E [exp(i {kiXt, +... + kiXt, + wMt,))] (2.54) 

({Hkj + ivistyyiHkjmy^-^ \ .3551 

IN ^'^\ {i){kj + w\6t)p {v{w\5t)y-'^ otherwise J ^ ' ' 

where in 1^2.55}) we replaced the delta functions in the integrand by exponentials 
eyi^{ikjXj) and exp(iu'j/) and we we factored out expressions depending on £,p. 
Example: 

(1) We take I — 1 and we obtain the joint characteristic function of the log- 
price and of the maximum of the log-price. We have: 

E[eM^{kXt + w^m = l(^(fe|<5t))'^ ( '^^^^ ^ J (2.56) 

The result 1^2.55]) will be used in in future work for pricing exotic options and 
for deriving variants of the Wiener-Hopf factorisation formula \21!/ for the operator 
stable Levy process in question. 



2.5. The option price 

An option on a financial asset is an agreement settled at time t to purchase (call) 
or to sell (put) the asset at some time T (maturity) in the future. For the sake of 
simplicity we consider European style options, ie such that can be only exercised at 
maturity (this means that boundary conditions are imposed on the option price at 
maturity t — T). Extending the analysis to American style options (to be exercised 
at any time ) is possible and can be done by considering European style options with 
m different maturities \14\ and performing the limit m —> oo. 

In order to minimise the risk we diversify the portfolio by dividing the money 
available in to Ns stocks St and Nc options C{St',t). This means that the portfolio 
V(t) reads: 

V{t)=NsSt+NcC{St;t) (2.57) 

where the coefficients Ns and Nc are customarily chosen \15'l as Ns = —dC /dSNc 
and in the following we take Nc = 1- 
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The investment strategy consists in replicating the portfolio, ie investing a copy 
of it into a safe central-bank account and let it earn a compound interest at a rate 
r independent oft. We analyse the distribution of deviations 

S)t := F(t + (5t) - e''(**) V(i) (2.58) 

between the compound interest e'"^**-''^ V^(i) — V(t) earned in the safe central-bank 
account and changes V{t + St) — V(t) of the price of the portfolio and we check 
if a self-financing strategy exists, ie if it is possible to choose C — ClSt^t) subject 
to a condition Ct = max{ST — K, 0) such that the average value of the portfolio 
deviation equals zero 

E[^i] = Q (2.59) 

(the deviations have no drift) and the variance of this deviation Var\I)i\ is minimal. 
Since we assumed 1^2.6]) that the waiting times between stock price changes are power 
law distributed with index /? we need to take the infinitesimal interest earned in time 
period [t,t + 5t] equal to r{St)^ in order to ensure that the deviation variable Di 
exists in the limit St ^ 0. 

Our approach differs here from the approach in financial mathematics where 
one assumes in the first place that the deviations St have no drift (there are no 
arbitrage opportunities on the market) and one imposes conditions (semi-martingale 
condition) on the high end of the jump intensity vj{x) ~ e~'^' in order to ensure this 
assumption. We cannot make such an assumption for operator stable Levy processes. 
The processes that we are using are sums of a large number of of independent, 
identically distributed jumps where the distribution of the latter diminish according 
to a power law and not to an exponential in the high end. The processes are invariant 
under auto- convolutions 1^2. 8\) (self- similar) and in our opinion it is much more 
appealing from the physical point of view to use them in financial modelling rather 
than to use the abstract notion of the metric space of Levy processes with right 
continuous and left bounded paths endowed with a Shorokhod topology if there is no 
experimental justification that financial time series in question are martingales. 

Our approach is not equivalent with theory in chapters 20.4-3 - 20.4-5 pages 
I4I6 - 1428 in }19f . That theory assumes that the log- characteristic function of 
the fluctuations that drive log St is analytic ( is power expandable ) and derives a 
partial differential equation (Fokker- Planck equation) for the price of the option as 
a function of the price of the stock and of time and solves this equation subject to 
initial conditions for the price of the stock at maturity. Operator stable processes are 
not power expandable. It might be possible to expand the log- characteristic function 
■uj{k) in a series of fractional powers of k using the fractional calculus ^221 and in 
particular the fractional Taylor expansion (formula (4-11) po-g^ 89 in ]2i^) \2'^ (^nd 



^We understand the price formation process as an accumulation of a large number of shocks that, 
once they exceed some threshold value, give rise to a change of the stock price (see discussion 
under Table □ 
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develop a fractional Fokker-Planck equation for the price of the option on the stock. 
In future work we will pursue investigations in this context in order to reconcile our 
results with that from flf^ . 

Equation 1^2. 1\) defines the price of a stock in terms of an exponential from a 
Levy process. Thus a small change of the stock price reads: 




[(d-5Lt)]-l) (2.60) 

("^ 



exp 

■5Lt 



E 

m— 1 



(2.61) 



where in 12.61\) we expanded the exponential in a Taylor series and retained all 
terms of the expansion since all these terms are proportional to (5t)^ 1^2. 2U\) for 
large values of 5Lt . 

From \2.6U\) and ^2.61]) we see that the price of the stock St is not a martingale. 
It is also not possible to construct a martingale from the price of the stock by sub- 
tracting the sum of means E [St+st — St\St] over t (the compensator) from St OS'- 
Indeed, the mean value of the left hand side of ^2. 6U\) conditioned on the history of 
the stock price up to time t is a sum of all possible positive powers of the fluctuation 
term ( a ■ SLt J and is, on the grounds of \2.19\) . infinite. It is therefore readily 
seen that, in the framework of our model, satisfying the condition \2. and thus 
solving the option pricing problem in its classical formulation may not be possible. 
We will therefore analyse the probability distribution of the deviation variable St 
and work out conditions under which it is possible to minimise the drift E [St] of 
the portfolio jg.57| ). 

We expand a change of the price of an option C{St', t) on the stock in a Taylor 
series in powers of the change of the stock price: 

Ct^s. - Ct = '-§^iStf + g(5,,, -St) + t (2-^2) 

n=2 

where the partial derivatives dC /dt an d^C /dS^ are assumed to be independent of 
time t and they will be determined later. 

From Y2.fi(]\) and \2.61]) we compute the nth powers of the changes of the stock 
price and we retain only those terms that are proportional to (St)^ . We have 

{St+st - StT = SI' (exp [(a • SLt)] - l)" forn>2 (2.63) 

/ ^ \ "2 ____ 

where we have neglected all terms [jiSt)^^ I a ■ 6Lt I on the grounds of \2.2(!l\) . We 
insert Y2. 621) and \2.6^) into \2.5t^) and we get the deviation random variable: 

S)t = [dtC + rStdsC - rCt) {Stf + ^ ^^^t" (exp [a ■ 6Lt) - 1) (2.64) 

Note that all terms on the right hand side of \2. ()4\l o'^s proportional to {St}^ . Now 
we derive the Fourier transform of the probability density function ofDt conditioned 
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on the value of the price of the stock at time t ( conditional deviation characteristic 
function). We have: 

X^,\sM:^E[e^'^^^\St]= (2.65) 

/°° / °° 1 f)nn \ 
dzv^^\z\5t)eyiY> ifc^ —J^^t (exp (z) - 1)" 42.66) 

/oo 

d\v^^'>{\\5t)m{k,X) ^ (2.67) 
-OO 

poo -» 

ey.Y>{ikDt{6tY) / d\w{{5tf&^ — {\))Tl{k,\) (2.68) 



where Dt := {dtC + rStdsC — rCt) and 

1 9"C 



m{k,X):=— dzexpUfc^ 

\ n=2 



5*4" (exp (z) - 1)" exp(-iAz) (2.69) 



anc? in 1^2. 61^} we made use of i2.15\) and of ^4.10}) . In further calculations we make 
use of the following proposition. 

Proposition 2.3. The integral 

1 1"°° 

m{k, A) := — / dz exp {ikP{e' - 1)} e"'^" (2.70) 

where P{z) is a real function such that the only root of P is P{0) — and P{oo) — 00 
tends, for small values of k, to a following delta function: 

Um,^,^{k, X)=5[x+ ^^^^^p'^_,l^ ^^^^ (2.71) 

where m > is a number that satisfies 1^5.11}) . P^^^^{z) is the inverse function to 
P{z) in some vicinity of z — and e > is a positive real number. The proof is 
given in Appendix B. 

Inserting \Z. 7j| ) into 1^2. we get the conditional deviation characteristic func- 
tion in the limit of small k. We have: 

X^^is, (k) = cM^kD.iStf)^ (^'^^)"^^M^ (log(P(-'^.0(f))) ^) ^'-^'^ 
where 

P{x) y ^^^r^^" (2.73) 

and we check by direct substitution that for d = 1 and E_ — the class of functions: 

°° 1 r)"r^ 

P^^)--=ll-l^Stx'' = ^(log(x))" (2.74) 

n=1 



February 7, 2008 14:52 WSPC/INSTRUCTION FILE OptionPriccCou- 
pledPriceTime 



Instructions for Typing Manuscripts (Paper's Title) 17 

with < a < ^ and T £ R and P{0) — ensures that the drift of deviations 



dk 



= (2.75) 



fe=0 

is zero. The equation \2. 15\j yields a differential equation for the option C = C{St, t) 
as a function of the price of the stock St and of time t viz 



m 



= dtC + rStdsC - rCt -Umk^Q ^ ' 



C2 \ / 1 p(-l,*) 

dtC + rStdsC + -^ds^C -rCt]-[ —limk^o- ^ ' 



2" 7 UmA' ^■^"(p(-i.»))'(f) 




Note that the first term enclosed by parentheses in \2. 7^ corresponds to the classi- 
cal Black& Scholes equation and the second term depends via \2. 74{ l on all partial 
derivatives ds"C. In future work we will analyse the uniqueness of the solution of 
the option pricing problem and we will derive conditions for the solution to have: 

Ct = C(5t, T) = max{ST - K, 0) (2.79) 

a given initial value at maturity T . 

3. Conclusions 

We have applied the technique of characteristic functions to the problem of pricing 
an option on a stock that is driven by operator stable fluctuations. We have shown 
that it is possible to construct an option on a stock as a function of the price of 
the stock in such a way that deviations of the portfolio are compensated by the 
compound interest earned on the portfolio. In future work we will price exotic options 
with different exercise times by computing the unconditional, joint characteristic 
function of the deviations of the portfolio related to values of stock price at these 
exercise times and by using initial conditions on the value of the portfolio at these 
exercise times. Testing of the theory in Monte Carlo simulations and fitting the 
theory to market data will also be pursued in future work. 

4. Appendix A 

We prove formula 1^2. 15\) for the conditional n-marginal probability function: 



i)^"-\k\6t)= e'''''i^^"-\z\6t)dz^ (4.1) 

J ^ expzfc (ct • SkYiJsUiST {^^t\^i) d{Slt) = (4.2) 
expifc (ct • Sk^ "(j) (sit\St^ d{Slt) = (4.3) 
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[dt) = I expi 



k (a ■ SltY vD{{5t)-^Mt)d{5lt) 



-f3TrE 



dXw{X) j 5lt exp(2 (k (a ■ Sk^ " - {{Sty'^M^'x) ■ Sk^ ) 



(St) 



-pTrE 



dXw{X) ( d{Slt)exp{i(t{Sky -l-Sh 

.lad \ 



(27r)° 



dXw{X)(Y{2Ti5{U)) / dw;e*'('""-"')e*('-'^)'" = 



J fI 27r J 

— oo 



(4.4) 
(4.5) 
(4.6) 
(4.7) 

(4.8) 

(4.9) 



In we conditioned on the fluctuation 5lt G K.'* and we integrated over 2 G M. 
In \4--S^ we assumed that the time St is small enough that at most one jump occurs 
during it and therefore we replaced the conditional fluctuation pdf by the condi- 
tional jump pdf and in {4-4}) we used \2. ?| ). In \4-^ we replaced the jump pdf w 
through an integral over X € K'^ from the Fourier transform w and we changed the 
order of integration integrating first over X and then over Sit ■ In | |^.6p we defined 
t /ccr", I :~ ^{{5ty^= X) we introduced an orthogonal operator ^ such that 

= (|(t|, 0, . . . , 0) and we substituted for ^6lt. In j^.Tj ) we performed the integra- 
tions over {5lt)2, ■ ■ • , {Slt)d o,nd in the remaining integral over {5lt)i we substituted 
for w = {Slt)i- In i4-'^ we eliminated the delta functions and transformed the d- 
dimensional integral over X into a one- dimensional integral over li . In ^4-9^ we 
defined a function (kernel) as : 



oo 



2Re 



(tll)'^ — Iw) 



" Jo '^^e 
2e*^r d^e"'^" cos ^cos(7r 



if n is odd 
if n is even 



(4.10) 



(4.11) 



and we evaluated it by rotating the range of integration in the complex plane by 
an angle 7r/(2n) (Wick rotation). In doing that we applied the Cauchy theorem 
to a contour composed of an interval w G [—R, R] in the real axis, two arches 
w = i?e** for (j> G [0,7r/(2n)] and G [tt, tt + n/{2n)] respectively and a line w = 
^gi7r/(2n) ^ g [—R,R]. The integrals over arches are bounded from above by 

R /J^^^^"^ (i(/)exp [— fii?" sin(n(/)) ± [i?sin((/))] for the arch in the upper f+J and the 
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lower (^) complex half-plane respectively and they tend to zero for R oo if n is 
even and positive. 

Now we compute a small ? expansion of the kernel /C^"-' [) for n even. We 
denote an '■= cos(7r ^"2""'"^ ), we substitute z := 6^" in j^- j j| ) and get: 

/CH(M):=;^/^rf-^/"-exp(-z)cos(|^.V«) = (4-12) 

e^-i-irf^V"' rdzz(^-)exp(-z)= (4.13) 



m=0 



?i7rm Vf/ ^ (277i)! ^ ' r- 



m— 



/n j^. we defined y (m7r/(2a„Z))" I, we expanded the cosine into a Maclaurin 
series and, owing to the fact that the period of the cosine is small for small t, the 
integrand oscillates rapidly and only the vicinity of the singularity of the integrand 
contributes meaningful to the integral, we truncated the integration at z — m7r/2 
for some m £ N. In i4.14)) we expressed the result through the truncated Gamma 
function. In j^. j.5| ) we defined 

C°° ( — 1)™ 7^2)71+1 \ 

S(2-+i)!nu(^+'^)j ^'■''^ 

we made use of a series expansion j^. j7| j of the truncated Gamma function. 

T{a,^)^^-e-^Y. f ^ ^ (4.17) 

We note that for t small enough the kernel satisfies: 

^ (4.18) 

for i = q^\q = (m7r)/(2a„gi/") • (t)^/" and for q-^ = q{m)-^ e N and q < 1 and 
thus we have 

/zme^o/C(") (t, l)c^s(l- [tf A (4.19) 

where £ := (m)/(g(m)"'^/"') is some number. It is readily seen from \4-l^ '^i^d j^. i j| ) 
that dl/C(") {i, [) = 1 and the same holds for the right hand side of ^4-^^ - 
Inserting into ^.Uj ) we obtain \2. 1 7| ) q.e.d.. 



February 7, 2008 14:52 WSPC/INSTRUCTION FILE OptionPriccCou- 
pledPriceTime 



20 Authors' Names 



5. Appendix B 

We calculate the integral \2.70^ for small values of k. We have: 



1 r°° 

A) := — / dz exp {ifcP(e" - 1)} e"*^^ = 

1 >^ r*+^ dw (p^-'^'Hf))' ™ / 

> / —-r^ — rT^¥^e™exp -lAloe 

^ I k [1 + P(-i).» (H)] ^ V ^ 



(5.1) 



27r 



1 ^ p^d^ (P(-^)--(f)) 



(-«Alog 



w 



1 

2^ 



E 



[P(-i)-^(f)] 



1+iX 



p(-i),i(0) z 



1+iA 



exp(-zAlog (p(-i)^* (H))) 



(5.3) 

(5.4) 

(5.5) 

(5.6) 

(5.7) 
(5.8) 



In \5.^) we substituted w = kP{e^ — 1), we denoted by Zi{y) := P(~^'''(y) the 
i^^ solution (inverse function) of the equation P(e^ — 1) — y and we divided the 
integration range over z into intervals [ai,ai+i] in which the function P(e^ — 1) 
is monotone and thus invertible. We also denoted ri := kP {e°'^ ~ !)• il.5.,Vj) we 
used the fact that P(oo) = oo and thus P("i)(oo) = oo and we approximated the 
integrand by its large w/k values. 

In {5.4)) we rotated the integration line by 7r/2 in the anti- clockwise direction 
making use of the fact that the integral over a quarter of a circle w — Pe**^ for 
(j) (z [0, 7r/2] is bounded from the above by J^^'^ d(f>e~^^™^'^^ and it vanishes for 
R — > oo. In i5. 5^} we used the fact that the only root of P'^^^^z) is at z — 
and thus we separated the integral into two parts the first one equal to ^ — 
over the integration range with all singularities excluded il = (IJ^ [r^, Ti+i] \ [0, e]) 
and the second one over [0, e]. Note that since Zimfc^oki+i =0 and the in- 
tegrand has no singularities in il the number S^/t ^ when k 0. The num- 
ber 9 counts intervals [r^jr^+i] such that G [ri,ri-|_i]. In J^5.5}) we substituted 
z = p(^i)^* (^) and in |5.6]) we performed the integral. Now we take m > 0, small 
k and Am ■— (— im^^) / (log(p(~^''^ (t)) '^^'^ evaluate the integral 1^5.8}) for this 
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value. 

mik, Am) + ^9^^^^ ^^k~ mexp {~m-') 6 (log(p(-i-') (f 

(log(P(-i,.)(^)) 

We fix k and let m — > 0. Then Am and we see that OT(/c,0) — > 0. Furthermore 
from 1^5.1]) we have dA*yi(fc, A) = 1 we conclude that 

Um,^onk, A) = 5 (^A + (5.10) 
and the number m can &e evaluated from 



|Am|^(fc,Am) = e-""^= i (5.11) 



q.e.d. 
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